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ABSTRACT 

We examine contributions of second order physical processes to the results 
of stellar evolution calculations which are amenable to direct observational test- 
ing. In the first paper in the series (Young, Mamajek, Arnett, & Liebert 2001) 
we established baseline results using only physics which are common to mod- 
ern stellar evolution codes. In the current paper we establish how much of the 
discrepancy between observations and baseline models is due to particular ele- 
ments of new physics in the areas of mixing, diffusion, equations of state, and 
opacities. We then consider the impact of the observational uncertainties on 
the maximum predictive accuracy achievable by a stellar evolution code. The 
sun is an optimal case because of the precise and abundant observations and 
the relative simplicity of the underlying stellar physics. The Standard Model is 
capable of matching the structure of the sun as determined by helioseismology 
and gross surface observables to better than a percent. Given an initial mass 
and surface composition within the observational errors, and no current observ- 
ables as additional constraints for which the models can be optimized, it is not 
possible to predict the sun's current state to better than ~ 7%. Convectively 
induced mixing in radiative regions, terrestrially calibrated by multidimensional 
numerical hydrodynamic simulations, dramatically improves the predictions for 
radii, luminosity, and apsidal motions of eclipsing binaries while simultaneously 
maintaining consistency with observed light element depletion and turnoff ages 
in young clusters (Young et al. 2003). Systematic errors in core size for models of 
massive binaries disappear with more complete mixing physics, and acceptable 
fits are achieved for all of the binaries without calibration of free parameters. The 
lack of accurate abundance determinations for binaries is now the main obstacle 
to improving stellar models using this type of test. 

Subject headings: stars: evolution - stars: fundamental parameters - hydrody- 
namics - convection - lithium 
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1. INTRODUCTION 

Stellar evolution has become a successful tool for elucidating the processes at work in 
individual stars. New instrumentation and a wealth of new data has resulted in increased 
emphasis in astronomy on the evolution of galaxies over cosmic history. Obviously the 
stellar content of a galaxy plays a central role in its evolution. In order to understand this 
process, we require theoretical stellar evolution to be predictive, in the sense of being able 
to accurately describe the contribution of luminosity, kinetic energy, and nucleosynthetic 
products from a star of a specific initial mass and composition at any and all points in its 
life. This process must be able to be carried out self consistently for stars from the hydrogen 
burning limit to the highest possible masses, so that stellar outcomes can be reliably linked 
to an initial mass function (IMF). This is not possible with schemes which are adjusted to 
match astronomical observations. Without an underlying physical theory, the calibration 
cannot be reliably extrapolated to regimes without extensive and independent observational 
data. Also, multiple physical effects can contribute in opposite or orthogonal senses to the 
star's structure. As an example, determinations of metallicity of binaries are often made 
by fitting model tracks of varying composition to observed points and finding the best fit. 
However, more mixing produces more luminous models for a given mass at a given effective 
temperature. Lower metallicity shifts tracks bluewards and slightly higher in luminosity. 
Either change could force a model to pass through the desired observed point. So a model 
with incomplete mixing physics and solar composition could be as good a formal fit to specific 
observations as a model with more complete physics and supersolar metallicity. 

The two primary areas which strongly affect the evolution and have uncertain physics are 
mixing and mass loss. The evolution is also sensitive to the opacity of the stellar material, 
but the opacities used in evolution codes are largely uniform, utilizing the OPAL values 
of Iglesias Sz Rogers (1996) for high temperatures and Alexander & Ferguson (1994) for 
low temperatures. The observational errors in determining stellar compositions are the 
major limitation on further testing contribution of opacities to stellar models. Even the 
metaUicity of the sun varies from z=0.019 (Grevesse & Sauval 1998) to z=0.015 (Lodders 
2003), depending upon the precise physical assumptions and dimensionality of the codes 
being used to fit the measured line profiles. Most other test cases, such as most double- 
lined eclipsing binaries, have no published metaUicity determinations using high resolution 
spectroscopy. Equations of state (EOS's) are not uniform across stellar evolution codes. 
While the effects of the EOS are perhaps more subtle, they can still be important, particularly 
for low mass stars and pre- and post-main sequence evolution. 

In Young, Mamajek, Arnett, & Liebert (2001), we presented baseline results from stellar 
models calculated using only physics common to current widely used stellar evolution codes. 
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These models were tested against a subset of double- lined eclipsing binaries (Andersen 1991; 
Ribas et al. 2000; Latham et al. 1996; Lastcnnct & Valls-Gabaud 2002; Hillenbrand & White 
2004). Young et al. (2003) discussed hydrodynamic mixing within the radiative regions 
of stars and presented several observational tests of the mechanism. This paper presents a 
reanalysis of the eclipsing binary sample and solar models, with more realistic mixing physics 
as well as additional minor improvements to the code. Section 2 summarizes the additional 
physics and improvements to the code. Solar models are examined in Section 3. The eclipsing 
binary sample is presented in Section 4. Section 5 contains discussion and conclusions. The 
implications for post-main sequence evolution will be presented in a subsequent paper. 



2. THE TYCHO CODE 

The TYCHO code is a ID stellar evolution and hydrodynamics code written in struc- 
tured FORTRAN?? with online graphics using PGPLOT. The code is as described in Young, 
Mamajek, Arnett, & Liebert (2001). We have made substantial additions and improvements. 
The code is now functional for stars from the hydrogen burning limit to arbitrarily high 
masses, and for metallicities of z = to the hmit of the OPAL opacity tables. 

The opacities used are from Iglesias & Rogers (1996) at high temperatures and Alexander 
& Ferguson (1994) for low temperatures. The OPAL tables have been extended to low 
entropies, and are formally adequate for calculating stellar interiors down to the hydrogen 
burning limit. In reality, a number of contributions, particularly from molecular species, are 
not included. Stellar models computed with these tables are reliable to perhaps 0.5 M©. 

TYCHO uses an adaptable set of reaction networks, which are constructed automatically 
from rate tables given a list of desired nuclei. In these calculations a 1?6 element network 
complete through the iron peak was used at T > 10^ K, and a 15 clement network for light 
element depletion at lower temperatures. Rates for the full network are from Rauschcr & 
Thielemann (2000). Caughlan & Fowler (1988) rates are used in the light element network. 

Mass loss capabihties of the code have been extended. At Tg// > ?.5 x lO^K, the 
theoretical approach of Kudritzki et al. (1989) is used. At lower Te// routines based upon 
the empirical prescription of Dupree & Reimers (1987) or Blocker (1995) are available. Low 
temperature mass loss was not important in any of the cases studied here, and the Reimers 
and Blocker algorithms converge in the limit of low luminosity. A treatment of radiatively 
driven mass loss in Wolf-Rayet stars based upon the work of Lamers & Nugis (2003) is also 
included in the code. It does not come into play for these models and will be discussed in a 
separate paper. 
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The equation of state has been expanded from the modified Timmes & Swesty (2000) 
EOS in Young, Mamajek, Arnett, & Liebert (2001) to include a more generaUzed treatment 
of the coulomb properties of the plasma. The formation and dissociation of molecular hy- 
drogen and its effect upon the equation of state are also included in a Helmholtz free energy 
formulation. The OPAL project has extended its EOS determinations to lower entropies. 
The improved TYCHO equation of state agrees with the OPAL EOS to better than 1% for 
most conditions. There remains a difficult region (10^^ ^ P ^ lO^g cm^ and T < 10^'^ K) 
in which the plasma is a strongly interacting coulomb system, and in which the difference 
exceed 4%. This region is relevant for low mass stars (M < O.5M0). 

TYCHO uses a modified version of Ledoux convection which avoids the problem of 
instantaneous mixing in convective regions with nuclear burning during short timesteps. 
Simple Eddington- Sweet rotational mixing (Tassoul 2000) is implemented in the code. Ro- 
tation was not included in the binary models as these systems are all relatively slow rotators 
for their mass range. Based on test models, rotation seems to be a perturbation of the order 
of the observed errors or less on the observed quantities we are examining. 

All models were run with an improved version of the inertial wave driven mixing de- 
scribed in Young et al. (2003). This treatment is primarily derived from Press (1981) and 
references therein. Press begins from two separate approximations to describe the propa- 
gation of internal waves in a stellar interior. The global behavior of the waves is described 
by the anelastic approximation (Ledoux 1974). Boussinesq, plane equations with thermal 
diffusion describe the damping behavior of linear internal waves (Lighthill 1978). Press 
(1981) combines these equations to define a set of globally anelastic but locally Boussinesq 
equations. The principal equation of motion is 
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N is the composition dependent Brunt- Vaisala frequency. 
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kn is the local horizontal wavenumber, 



kn = 



(4) 



a is the thermal diffusivity, 
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and uj is the characteristic frequency of the driving, in this case convective modes of 
low frequency. The vertical fluid velocity uy (for moderate damping) is determined from a 
WKB solution of equation 1, and the vertical wavenumber ky is defined as 



ky = kn 
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The quality factor Q of the propagating waves is 
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Waves generated by the convective zone are in most circumstances non-linear within 
one linear damping distance of the convective zone. Mixing in this region is efficient, but 
it is a small fraction of the star. We are also interested in the impact of the waves outside 
this zone. Damping of the waves in the radiative region generates vorticity as described in 
Young et al. (2003). Press (1981) derives a mass diffusivity coefficient gm, which we adopt, 
as. 
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where e is the dimensionless nonlinearity parameter for internal waves 
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The efficiency of the wave propagation, and therefore the mixing, is inherently sensitive 
to composition gradients through N and, indirectly, through the opacity in a. The mixing 
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occurs roughly on a thermal timescale. It is slow relative to convective turnover, but it can 
be significant on evolutionary timescales. In practice, this drives mixing in what are typically 
referred to as "semiconvective" regions, where — 1) > 0. When this quantity becomes 
imaginary, the wave undergoes complete internal reflection, and there is an exponential falloff 
of the wave flux beyond this boundary. In terms of composition gradients, there is a region 
outside the convective zone which is well mixed, with a further sharply decreasing degree of 
mixing beyond. This treatment removes the one variable parameter from Young et al. (2003), 
the characteristic length for dissipation, by treating the wave physics more completely. Our 
results are little changed from the earlier formulation. This is unsurprising, as we constrained 
the earlier treatment with numerical simulations in which the wave behavior was readily 
apparent. 

Gravitational settling and differential diffusion of nuclear species according to Thoul, 
Bahcall, & Loeb (1994) is also included. The Thoul, Bahcall, & Loeb (1994) treatment 
of diffusion is gcncralizable to an arbitrary number of nuclear species, though that work 
examines only H, He, O, and Fe. We calculate diffusion coefficients separately for the species 
important to the OPAL "type 2" opacity tables (H, He, C, and O) and on average for iron 
peak and Ne like elements. Michaud et al. (2004) examine the effect of settling on the 
approximately solar age and metallicity clusters M67 and NGC 188, using 19 elements. In 
light of these results our intermediate simplification appears adequate for stellar structure 
calculations. 

Numerous minor improvements have been made which improve convergence and stability 
of the code, and allow it to perform adequately at the small timesteps typical of neutrino- 
cooling dominated burning stages as well as the slow hydrogen burning stages. The code 
is publicly available and open source. The current version (TYCHO-7.0) is being made avail- 
able, along with an extensive suite of analysis tools, at http:/ /pegasus. as. arizona.edu/~darnett 

3. SOLAR MODELS 

As the best observed star in the sky, the Sun is an obligatory test case for any com- 
prehensive stellar evolution code. The hehoseismological measurements of sound speed and 
depth of the convective zone give us an insight into the interior structure not available for 
any other star. In this section we test solar models from TYCHO, but with a somewhat novel 
aim. We hope that TYCHO will function as a predictive tool for building stellar populations. 
Therefore, instead of finding a best fit to objects as they are observed now, with variable 
initial conditions, the code must be able to predict a unique (and accurate) path through 
stellar parameter space over time for a particular initial mass and composition. Conversely, 
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we would also like to connect any given observed star to a unique initial condition. We would 
wish to do this for the complete range of stellar masses. As such we are more interested in 
the comparison of our models with the sun, assuming only an initial solar mass, composition, 
and our best treatment of the physics, than in how precisely we can fit the sun by optimizing 
our models. A 1 Mq star on the main sequence is probably the easiest type of star to model, 
being relatively insensitive to the effects of mixing and mass loss. Solar models give us an 
estimate of the minimum uncertainty in our predictions of stellar parameters. 

We examine four models, sO, si, s2, and s3 which differ in the completeness of mixing 
physics included. Model sO includes gravitational settling and diffusion (Thoul, Bahcall, & 
Loeb 1994) and inertial wave-driven mixing (Young et al. 2003). Model si includes only 
wave-driven mixing and s2 only gravitational settling and diffusion. Model s3 uses only 
Ledoux convection and ignores other mixing physics. Eddington-Sweet mixing is disabled, 
as it is a poor description of the true angular momentum distribution in the sun. We also 
calculate one model (10) with physics identical to sO, but with Lodders (2003) values for solar 
abundances. 

There is one glaring limitation, a free parameter which must of necessity remain in this 
1-D code. We choose a mixing length parameter of a = 2.1, where a. is the ratio of the mixing 

length to the pressure scale height. This is in the same range as values deduced from solar 
standard models [a = 2.05) (Basu, Pinsonneault, & Bahcall 2000) and multidimensional 
simulations of the solar convective zone with hydrodynamics and radiative transfer {a = 
2. 13) (Robinson ct al. 2004). Smaller values of the mixing length parameter result in larger 
radii for the 1-D models. The mixing length parameter is a shorthand for more complex 
physics, and is not guaranteed to be constant over stellar mass or evolutionary state, so 
these constraints apply only to sun-like stars. 

Table 1 gives values for gross obscrvablcs {R,Teff,L), the rms difference in predicted 
and observed sound speed, depth of convection zone, photospheric He and Li values, and 
central temperature for each model. (Note that the observed Li abundance in the sun is 
much lower than the meteoritic value. The meteoritic value is assumed to more accurately 
represent the abundance in the proto-solar nebula, and as such is used as the starting value 
for the models.) We have not performed an inversion of the helioseismological data through 
our model to obtain expected sound speeds for our models. The values to which we compared 
are those calculated by Bahcall, Pinsonneault, & Basu (2001). We have roughly twice as 
many Lagrangian mass points as the online tables for the standard models, so we perform 
an interpolation in our sound speeds to match the standard model tables. Performing this 
direct comparison we find rms errors in sound speed of ~ 0.5% for our best models. Most of 
this discrepancy can be attributed to our error in the solar radius. (Bahcall, Pinsonneault, 
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& Basu (2001) find a 0.15% rms error for a model with a 0.04% difference in radius from 
their standard value.) 

Table 2 gives neutrino fluxes for the models and a selection of models from Bahcall 
& Pinsonneault (2004). The Bahcall & Pinsonneault (2004) model BP04 is the current 
standard solar model. BP04+ includes improvements in nuclear physics, solar equation of 
state, nuclear physics and solar composition. The entry "Comp" shows the effects of variant 
composition alone, and N^^ shows the effects of a change in the N^^(p, 7)0^^ rate alone. 

The values in Table 1 illustrate some of the subtleties involved in distinguishing between 
models. If we accept a constraint on the mixing length from simulations or helioseismology, 
all variants of the model predict gross observables to within 3%. The models with more 
complete mixing physics show a slightly better agreement, but the variation is less than the 
uncertainty in the exact nature of the convection. The minimum uncertainty in our pre- 
dictions must be take to be larger than 3%, because the error is dominated by a fictitious 
parameter. Varying the mixing length by 0.1 results is roughly a 1% change in the radius. 
Simulations of red giant atmospheres (Asida 2000) and observations of pre-main sequence 
(pre-MS) binaries (Hillenbrand & White 2004; Stassun et al. 2004) indicate that stars with 
low surface gravities and larger convective cell sizes and/or Mach numbers and turbulent 
pressures have different convective physics than main sequence stars of the same luminosity. 
In a ID description of convection, this manifests as a change in the mixing length to values 
of roughly 1.5. Without constraints on the nature of convection the minimum predictive 
uncertainty is roughly 7% for a IMq star of solar age. Varying the abundance from Grevesse 
& Sauval (1998) to Lodders (2003) introduces a further uncertainty of ~5% in surface ob- 
servables, and a much larger variance for hehoseismological and abundance tests. Model 10 is 
a poor fit primarily because the opacity at the base of the convective zone is too low. Press 
& Rybicki (1981) find that hydraulic enhancement of the photon diffusivity by wave motion 



Table 1. Characteristics of Solar Models with Varied Mixing Physics 



Model 


R/Rq 




log L/L0 


rms Scs (%) 


Rconv 1 R 




log e{J Li) 


Tc 




sO 


0.993 


3.765 


3.80 X 10-3 


0.56 


0.714 


0.242 


1.14 


1.560 X 


10^ 


si 


0.972 


3.770 


7.65 X 10-3 


0.90 


0.729 


0.279 


1.78 


1.564 X 


10^ 


s2 


0.985 


3.767 


5.93 X 10-3 


0.50 


0.718 


0.240 


3.05 


1.562 X 


10^ 


s3 


0.971 


3.770 


6.85 X 10-3 


0.97 


0.729 


0.279 


3.11 


1.563 X 


10^ 


10 


0.944 


3.777 


9.21 X 10-3 


2.24 


0.672 


0.252 


4.89 


1.558 X 


107 


standard modeP 


1.0 


3.762 


1.0 


0.10 


0.714 


0.244 


1.1 


1.569 X 


10^ 



^Solar values from standard solaj: model of Bahcall, Pinsonneault, & Basu (2001) except log e(J Li) from Boothroyd & 

Sackmann (2003) 
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acts in a radiatively stable stratified region to increase the effective opacity, so the Lodders 
abundances may not be irreconcilable with solar observations. This may also bear upon the 
perennial problem of sound speed discrepancies between models and helioseismology imme- 
diately below the solar convection zone. The implications of this process will be discussed 
in a subsequent paper. 

Using helioseismology and detailed chemical abundances, we can begin to discriminate 
between models. Unsurprisingly, model s3, with mixing limited to Ledoux convection, is 
ruled out immediately because its He and Li abundances are much higher than observed, 
essentially unchanged from their primordial values. Model si, with radiative region mixing 
but no heavy element diffusion, is also eliminated by the size of the convective zone and 
surface helium abundance. Clearly gravitational settling and diffusion are necessary to fit 
the observed sun. The only observable difference between the remaining diffusion only and 
more realistic mixing models lies in the predicted photospheric Li abundance. This is exactly 
what is to be expected, since helioseismology tells us that, while mixing must be present in 
the radiative regions to account for observed abundances, it cannot have a large structural 
effect. 

Michaud et al. (2004) confirm this result for the solar age and metallicity clusters M67 
and NGC 188. They find that models with little or no "overshooting" are consistent with the 
observed color-magnitude diagrams of the clusters. Our theory of mixing naturally predicts 
little structural effect on solar type stars and only a small increase in core size for stars 
which start with small convective cores. We do see a significant effect on the sun during the 
pre-MS, when the transient convective core established during partial CN burning is at its 
largest. 

With no mixing (save settling) outside the convection zone, model s2 greatly under- 
predicts the depletion of Li at the solar photosphere. The model with complete mixing gives 
an abundance much closer to the observed value of log t(J Li) = 1.1 ± 0.1 (Boothroyd & 
Sackmann 2003), though this, too is sensitive to the mixing length at the factor of 2 to 
4 level. The role of rotation coupled to oscillations in driving mixing has been discussed 
extensively by many investigators ((i.e. Chaboyer, Demarque, & Pinsonneault 1995; Pinson- 
neault et al. 2002)). The work of Charbonnel & Talon (1999); Talon, Kumar, & Zahn (2002) 
suggests that mixing is damped in rotating stars on the red side of the Li dip, corresponding 
to early G stars. If the pre-MS sun was a slow rotator, it may be a limiting case where 
angular momentum transport produces a minimal modification in the stellar g-mode oscil- 
lation spectrum and mixing is at a maximum. This provides a possible explanation for the 
strong depletion of Li in the sun relative to field G stars, and for the wide observed range of 
depletions, but the problem is beyond the scope of this paper. 
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The values for neutrino fluxes in Table 2 all fall within the range of variation found 
by Bahcall & Pinsonneault (2004) for variant models with improved physics. A selection of 
the Bahcall & Pinsonneault (2004) models illustrating the range of variation between the 
models are given in Table 2. The neutrinos do not provide a constraint on the models at 
this level, but do confirm that none of the physics included in the models is in confiict with 
the observations. 

The solar models highlight some of the problems in assessing the predictive power of 
stellar evolution codes. The models presented here would be indistinguishable for a G2 star 
outside the solar system. The errors in the gross obscrvablcs could be compensated for by 
a change in the mixing length without including the necessary physics of He and heavy 
element diffusion and non-convective mixing. Helioseismology would not be available to 
falsify such a model based upon convective zone depth or sound speeds. Stars with abundance 
determinations from high resolution spectroscopy are relatively few and far between. One 
may argue whether it is then important to have complete physics. Figure 1 illustrates the 
potential traps that lie in validating code with a narrow selection of observations. At the 
solar age models sO and s2 (both with diffusion) are nearly identical, as are si and s3 (both 
without). The effects of diffusion clearly dominate in determining the sun's position in the 
HR diagram. On the pre-MS during partial CN burning in the transient convective core, the 
case is very different. Models sO and si (inertial wave-driven mixing) are nearly identical, as 
are s2 and s3 (convective mixing only). The shape of the pre-MS is determined primarily by 
the change in convective core size resulting from including more complete mixing physics. 
Diffusion has had insufficient time to make much difference. In short, the evolutionary 
history is not unique. A model which fits the present day sun perfectly may be substantially 
inaccurate for other evolutionary stages (or equivalently mass ranges or compositions) where 
different physics come into play. Calibrations based upon any one type of data set should 
not be extended into other regimes unless based upon a valid physical theory. 



4. ECLIPSING BINARIES 

The double-lined eclipsing spectroscopic binaries in Andersen (1991) provide us with a 
sample of stars from IMq < M < 23Mq with precisely determined masses and radii. A 
subset of these stars also have measured apsidal motions, which provide some information 
on the interior density profiles and core sizes of the stars. Wc use the same sample as Young, 
Maniajek, Arnett, & Liebert (2001) so that a direct comparison of the same code with and 
without complete descriptions of mixing physics outside of convective regions can be made. 
We use the most recent values available for observed quantities (observed quantities and 
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Fig. 1. — Evolutionary tracks for 1 Mq stars with four variations on mixing physics. Model 
sO includes gravitational settling/heavy element diffusion and inertial wave-driven mixing 
in radiative regions. Model si contains only wave mixing. Model s2 includes only set- 
tling/diffusion. Model s3 has no mixing outside of convective regions. Models with hydro- 
dynamic mixing in radiative regions (sO and si) are indistinguishable on the pre-MS, when 
the shape of the track is controlled primarily by the size of the small convective core that 
exists during partial CN burning. Diffusion has insufficient time to work on the pre-MS and 
does not affect the tracks. On the main sequence, after the convective core disappears, the 
shape of the track is determined by the presence (sO & s2) or absence (si & s3) of diffusion 
and settling. 
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references can be found in Table 3.) 

As in Young, Mamajek, Arnett, & Liebert (2001), we calculate a x^-like quantity for 
each binary pair, defined by 

= ((logL(m^,t)-logL^)K(A))2 
+ {{\ogL{mB,t)-\ogLB)/aL{B)f 
+ {{log R{mA^)- log RA)/(yR{A)f 

+ {{\ogR{mB,t)-\ogRB)/aR{B))\ (10) 

where A and B denote the primary and the secondary star, respectively. Here La and Ra 
are the observationally determined luminosity and radius of the primary, with a la and ura 
being the observational errors in log La and in log Ra- We convert the observational data for 
the radii to logarithmic form for consistency. Correspondingly, L{mA, t) and R{mA, t) are the 
luminosity and radius of the model. This was evaluated by computing two evolutionary 
sequences, one for a star of mass tua and one for m^. A was calculated at consistent 
times through the entire sequence (t^ = to a fraction of a time step, which was a relative 
error of a few percent at worst). The smallest value determined which pair of models 
was optimum for that binary. Note that if the trajectories of both A and B graze the error 
boxes at the same time, x^ ~ 4. We use radius instead of effective temperature in our 
fitting algorithm because the more precise values for R make the x^ more discriminating. 
These error parameters along with the corresponding uncertainties from the observations are 
presented in Table 4. 

As before, we note that a x^ statistic assumes that the observational errors have a 
Gaussian distribution about the mean (Press et al., 1992). This is not necessarily true, 

as the systematic shifts in measured quantities due to new analyses can be much larger 
than the formal error bars (Ribas et al. 2000; Stickland, Koch & PfcifTcr 1992; Stickland, 
Lloyd, & Corcoran 1994). Also, the quoted luminosity depends on the radius and effective 
temperature, and is thus not entirely independent. These systematic errors are the true 
limit to our power to discriminate between models, and emphasize the need for independent 
observational tests and numerical simulations to identify relevant physics. 

Probably the greatest observational limitation we face is the lack of abundance deter- 
minations for these stars. The only binary in our original sample with a spectroscopic 
abundance determination is UX Men (2;=0. 019) (Andersen, Clausen, Magain 1989). A few 
other systems have some sort of metaUicity indication in the literature. Ribas et al. (2000) 
derive a metaUicity of 2;=0.013 from fits to evolutionary tracks for the ever troublesome ( 
Phe. Synthetic BaSeL photometry of VV Pyx suggests a metaUicity of z <0.007, but the fits 
are not good (Lastennet et al. 1999). Latham et al. (1996) and Lastennet & Valls-Gabaud 
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(2002) argue for a metal content similar to the Hyades in DM Vir (^=0.23). MY Cyg A&B 
and AI Hya A are all peculiar metal line stars. We (somewhat arbitrarily) also assign these 
systems a Hyades composition. Al Hya has a measured z = 0.07, but this is probably a 
surface enhancement and does not reflect the global composition of the star (Ribas et al. 
2000). Other systems either have no metallicity determinations or are sufficiently near solar 
composition that models of solar composition fall within the observational errors. 

4.1. Global Properties of the Errors 

The values for each binary pair with and without complete mixing physics are plotted 
in Figure 2. An arrow shows the shift in from baseline models to models with the full 
suite of physics. Fifteen of the binaries have excellent fits (x^ < 4). Only three systems arc 
marginal {( Phe, = 4.10; AI Hya, = 4.12; EK Cep, = 5.97). These systems will be 
discussed individually later. In all cases where the previous fits were marginal to poor, the 
improved. All massive binaries (M^, Mb > 4:Mq) with good fits also improved. Results 
were mixed for lower mass stars with good fits. In both of the latter groups, both complete 
and incomplete models fall within the observational errors, so there is little to distinguish 
between models for individual binaries. The threshold for rejection of a model with z/ = 4 
degrees of freedom with an a = 5% chance of rejecting a true hypothesis is = 9.488, so 
even the marginal fits do not give us a strong hold on remaining errors in our models. We 
must examine systematic discrepancies or wait for tighter observational error bars. 

We show values for our best estimate of the metallicity. It is difficult to make a 
global statement about how metallicity affects the results of the tests, because it varies 
depending upon the characteristics of the binary. For the most massive binaries, a change 
in composition from solar [z = 0.0189) to z = 0.004 (roughly that of the Large Magellanic 
Cloud) results in a shift in the HR diagram roughly comparable to the observational errors. 
A similar shift in abundances for the lowest mass stars can shift the models well outside the 
region of good fit, primarily because of the smaller size of the observational errors rather 
than an intrinsic physical sensitivity. The converse is not necessarily true, however. AI Hya 
is the most extreme case. Using a solar composition model with complete mixing physics 
results in a = 11-2, because the model is too luminous at the correct Tg//. The converse 
situation does not hold. A solar composition model with only standard physics has x^ > 20. 
A supcrsolar composition docs not produce a significant improvement because the main 
sequence never travels far enough to the red for either case. Only increased mixing solves 
this problem. In general, it is safe to say that a change in metaUicity of a factor of 2-4 is 
sufficient to move a binary model out of good agreement with observations. A similar change 
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may or may not make poor model fit well. In order to pin down the stellar physics we must 
either examine systematic trends or find better abundance determinations. 

Figure 3 shows goodness of fit vectors with the HR diagram for all stars in the sample, 
with incomplete models on the left and complete models on the right. The observed points 
with error bars are plotted with an arrow indicating distance and direction to the best fit 
model point. Wc can now begin to discriminate between models even for formally excellent 
fits. The most striking feature of the figure is the behavior of the massive stars (M > 4Mq). 
The incomplete models are systematically undcrluminous. This suggests three possibilities: 
(1) the massive stars arc all low metallicity; (2) the observational luminosity and/or mass de- 
terminations are systematically low; (3) the stars have larger convective cores than standard 
models predict. Option (1) is unlikely for nearby massive stars with ages of less than 10^ 
years, but cannot be absolutely ruled out without spectroscopic abundance determinations. 
Option 2 is possible, but again unlikely for a sample of 6 widely separated binaries. (Both 
EM Car and CW Cep have had their masses revised downward by Stickland, Koch & Pfciffer 
(1992).) Option 3 seems the most likely, and is consistent with other evidence of mixing in 
stars beyond the standard model. Indeed, the trend virtually disappears when realistic mix- 
ing is included. It may be argued that the standard model is hydrodynamically inconsistent, 
as indicated by the instantaneous deceleration required at convective boundaries, and by 
detailed analysis of realistic high resolution 3-D hydrodynamic simulations, problems which 
our mixing algorithm addresses. 

The salubrious effect of realistic mixing is confirmed by the apsidal motion tests. Models 
with incomplete mixing have systematically higher predicted apsidal motions in the four most 
massive systems. This indicates that the models are not sufficiently centrally condensed. (See 
Young, Mamajek, Arnett, & Liebert (2001) for a complete discussion of our previous results 
and methodology.) The current models with realistic mixing physics have larger convective 
cores and are therefore more centrally condensed. Table 5 summarizes the apsidal motion 
results for our binaries. Figure 4 shows the dimensionless rate of apsidal motion, {P/U)cl = 
{P/U)oBS ~ {P/U)gr, which would be due to classical apsidal motion, plotted versus log 
of half the total binary mass (where CL denotes classical and GR general relativistic parts, 
and OBS the observed motion). P is the orbital period and U the apsidal period. The 
observational data (corrected for general relativity) are shown as diamonds, with vertical 
error bars. All of the massive star models now fall within the error bars for the measured 
apsidal motion with the exception of QX Car, which differs by roughly two and a half sigma 
(insofar as sigma is a meaningful expression of these errors) . QX Car does not differ from 
the measured point by a larger absolute amount than the other binaries; it simply has much 
tighter error bars. Either the observational uncertainties are underestimated or, equally 
likely, the models are still missing some physics. We have reanalyzed all of the binaries from 
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Fig. 2. — values for optimum models of binaries, versus mean mass of the binary. Arrows 
are drawn from Young, Mamajek, Arnett, & Liebert (2001) values to values for current 
models with more complete physics. The improvement is dramatic. The vertical line at 
— 4: indicates a fit in which the models are just within the observational errors for both 
stars. 
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Fig. 3. — Goodness of fit vectors for each member of binary sample, with observational error 
bars. Young, Mamajek, Arnett, & Liebert (2001) results are on the top; present results on 
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EM Car to IQ Per. The lower mass binaries all have quoted apsidal motions smaller than the 
predicted general relativistic term. We find it more likely that there are errors in measuring 
an apsidal motion with periods of centuries than in weak-field general relativity. 

4.2. Individual Systems of Interest 

4.2.1. C Phe 

This system is perennially troublesome to stellar modelers (i.e. Ribas et al. 2000). It is 
difficult to fit both components with the same metallicity. We adopt z = 0.013, following 
Ribas et al. (2000) and achieve a marginal fit (x^ = 4.026). The secondary star is more 
luminous than the models when a good fit is achieved for the primary. If the observations 
are correctly interpreted, then the secondary model requires either a lower metallicity or 
enhanced mixing. Figure 5 shows the observed points and model tracks for ^ Phe. 

4.2.2. Al Hya 

Al Hya is identified as a peculiar metal line star (spectral class F2m). Its metallicity is 
measured as z = 0.07 (Ribas et al. 2000), but this is probably due to a surface enhancement. 
Still, the stars are probably metal rich relative to solar. Without a precise determination 
of the interior composition we choose to use a Hyades composition {z — 0.023) as being 
in the reasonable range of nearby metal rich compositions. In keeping with our effort to 
test the predictability of our code, we do not try to optimize the fit by further varying 
the composition. This system is particularly interesting in that increased metallicity alone 
cannot reconcile tracks with only convective mixing with the observations. The primary of 
the system lies farther rcdward in the HR diagram than the terminal age main sequence 
(Young, Mamajek, Arnett, & Liebert 2001). The lifetime on the Hertzsprung gap for a 2 
Mq star is short. It is possible to catch a star in that stage, but unlikely, as Tgap/rMs < 2%. 
More realistic mixing gives a larger convective core, extending the track redward so that a 
fit on the main sequence is easily achievable. Figure 6 shows the observed points and model 
tracks for Al Hya. 



Table 2. Neutrino Fluxes for Solar Models with Varied Mixing Physics 



Model 


PP 


pep 


hep 




B8 


Ni3 








IQlOcm^s-l 


lO^cm^s-i 


lO^cm^s-i 


lO^cm^s-l 


lO^cm^s-i 


lO^cm^s-l 


lO^cm^s-l 


lO^cm^s-l 


sO 


5.95 


1.42 


7.91 


4.83 


5.51 


4.08 


3.49 


4.59 


si 


5.97 


1.42 


7.84 


4.86 


5.59 


4.12 


3.53 


4.65 


s2 


5.96 


1.42 


7.84 


4.81 


5.44 


4.04 


3.45 


4.54 


s3 


5.97 


1.42 


7.85 


4.84 


5.51 


4.09 


3.49 


4.58 


BP04* 


5.94 


1.40 


7.88 


4.86 


5.79 


5.71 


5.03 


5.91 


BP04+ 


5.99 


1.42 


8.04 


4.65 


5.02 


4.06 


3.54 


3.97 


Comp 


6.00 


1.42 


9.44 


4.56 


4.62 


3.88 


3.36 


3.77 


Ni4 


5.98 


1.42 


7.93 


4.86 


5.74 


3.23 


2.54 


5.85 



^Neutrino fluxes from standard model and models with various improvements in physics from Bahcall & Pinsonneault (2004) 
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Fig. 4. — Classical apsidal motion versus mean mass, for our binaries with measured apsidal 
motion. {P/U)cl = {P/U)obs ~ {P/U)gr is assumed. Diamonds are observed values, 
squares are model fits. Note that the symbol for QX Car is larger than error bars, and QX 
Car is not a good fit. 
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Fig. 5. — Evolutionary tracks {z — 0.013) and observed points with errors for ( Phe. Note 
that the secondary star is more luminous than the models when the primary is well fit. 
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Fig. 6. — Evolutionary tracks {z — 0.023) and observed points with errors for AI Hya. 
Realistic mixing allows an acceptable fit to this previously difficult pair. 



4.2.3. EK Cep 



EK Cep and RS Cha are both pre-MS systems (Popper 1987; Mamajek, Lawson, & 
Feigelson 2000). The fit to RS Cha is formally a good one (x^ = 2.776), and we do not 
attempt to optimize within the observational errors. EK Cep, however, achieves only a 
marginal fit (x^ = 5.973), because the radius of the secondary star is larger that that of the 
models during the first rise of the CN burning bump. (Other early pre-MS models show 
similar behavior, so we suspect this is systematic. A larger sample will be discussed in a 
forthcoming paper.) This is a robust behavior, in the sense that most things wc coTild do 
to the models do not push them in the right direction. The non-convective mixing physics 
does not have a substantial effect, and increased metallicity would change the luminosity 
too much to result in a good fit. We can only find good agreement by reducing the mixing 
length parameter to a = 1.6. This suggests a change in the nature of the convection, but 
since a does not represent a physical entity, it does not tell us what that change is. We may 
speculate that since the star is trying to transport an amount of energy to a surface with a 
larger radius than a main sequence star of similar luminosity, the convective Mach numbers 
must be higher. This contributes a proportionally greater term to the stress tensor than 
main sequence convection and manifests as a radially directed pressure term, which would 
result in a larger radius for hydrostatic equilibrium. Besides this term, there are plasma 
effects, non-hydrogenic molecular contributions to the EOS, molecular and grain opacities, 
and subtleties of atmosphere models which must be taken into account which may contribute 
to the resolution of the problem. In short, we can identify a deficiency in our physics, and 
probably localize it to the physics of convection, but we do not have good predictive accuracy 
in this evolutionary stage. We quote the for our usual value of a — 2.1, and not the 
improved fit for a — 1.6, since this is not a predictable change. We plan 3-D simulations 
of convection in pre-MS and MS stars, which wc hope will characterize the difference in the 
convection in a physical way. Figure 7 shows the observed points and model tracks for EK 
Cep. 



4.2.4. TZ For 

TZ For was not in our original binary sample, but is a sufficiently interesting system 
that we examine it briefly here. The secondary is a subgiant in the Hertzsprung gap (Pols et 
al. 1997) with a spectroscopically determined metallicity oi z — 0.024 ± 0.007. Lastennet & 
Valls-Gabaud (2002) attempt to fit the secondary with several stellar evolution codes, but 
are unsuccessful without changing the mass of the model by 5a or using a composition not 
in agreement with the observations. Figure 8 shows the observed points and model tracks for 
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Fig. 7. — Evolutionary tracks {z — 0.019) and observed points with errors for EK Cep. The 
models for the secondary star were run with a reduced mixing length parameter of a = 1.6 
in order to increase the model radius to that of the observations. This indicates a change in 
the nature of the convection. 
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TZ For. 

Changing the size of the convective core on the main sequence necessarily changes the 
path the star takes across the Hertzsprung gap. We find that with reahstic mixing, our 
models match the hotter component of TZ For reasonably well, and the cooler component 
exceptionally well. The models for the subgiant are slightly overluminous. We find a, — 
0.77 for the binary at t — 5.3 x 10^ yr. Virtually all of this discrepancy comes from the 
subgiant. 

5. CONCLUSIONS 

In this paper we test the predictive power of the TYCHO stellar evolution code against a 
set of classical observational tests. With an improved version of the realistic mixing physics 
presented in Young et al. (2003), we find excellent agreement with solar models and the 
sample of double-lined eclipsing binaries from Young, Mamajek, Arnctt, & Licbcrt (2001). 
By avoiding optimization of our models with composition changes or parameterized extra 
mixing, we also identify several issues which are important to future development of stellar 
modeling. 

From the solar models we find that our predictive accuracy for the sort of gross observ- 
ables {L, R,Teff) we measure for other stars is limited to of order 5-10% by (1) inadequacy 
in our description of convection, manifested by an uncertainty in the fictitious mixing length 
parameter, and (2) by uncertainties in abundances. If the nature of the convection is fixed 
by numerical simulations of full 3-D convection, the uncertainty is reduced to that arising 
from the abundance determinations. The good agreement of the neutrino fiuxes with those 
of the standard model indicate that the infiuence of the mixing length description is an 
"atmospheric" effect. The rest of our (more or less parameter- free) physics provides a good 
description of the interior of the sun. We wish to emphasize that the quoted uncertainty 
is not a precise figure, but rather our best estimate, accounting for various known errors 
and incomplete physical descriptions. The sun is a well constrained system, and has un- 
complicated physics compared to most other stars. In this light our measure of uncertainty 
should be taken as a best case for stellar evolution, with the caveat that errors can become 
substantially worse as we examine more complex, less well observed systems. 

One of the most striking features of the solar models underlines a fundamental problem 
of stellar evolution. The two models that match observed solar quantities best have virtually 
identical tracks on the main sequence, which are shaped primarily by the inclusion of gravi- 
tational settling and diffusion of heavy elements. The models diverge significantly, however. 
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Fig. 8. — Evolutionary tracks {z — 0.024) and observed points with errors for TZ For. The 
models are slightly over luminous with respect to the observed hot component, but much 
improved over standard models. 
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on the pre-MS, where the influence of hydrodynamic mixing in radiative regions dominates 
the evolutionary pathway while the transient convective core is at its maximum extent. (Our 
theoretical treatment naturally predicts a smaller effect on the HR diagram for stars with 
smaller convective cores. This is consistent with the results of Michaud et al. (2004), which 
find minimal overshooting for the stars near the convective/radiative boundary in M67, a 
result which contradicts simple parameterized overshooting.) In at least some cases, a good 
fit to observations can be achieved without including physics which may be very important 
to the overall evolution. This may be adequate for describing the state of an individual star, 
but presents a serious problem for characterizing the behavior of a star or population over 
time. 

The eclipsing binaries provide a test of our physics, particularly the more complete mix- 
ing, over a wide range of stellar masses. The systematic problems with massive star models, 
which were identified in Young, Mamajek, Arnett, & Liebert (2001), are amehorated by the 
new treatment. The models are no longer under luminous, and the central condensations 
as measured by apsidal motions are no longer too small. Both of these improvements arise 
from larger convective core sizes resulting from the improved mixing. Simultaneously, the 
fits for almost all of the lower mass stars improve as well. The one case in which the error 
formally increases varies within the observational errors. All of the poor or marginal fits in 
Young, Mamajek, Arnett, & Liebert (2001) improve dramatically. Some of the improvement 
in these lower mass models arises from the use of non-solar abundances in a few cases, but 
the composition alone cannot account for all of the error in the earlier models. AI Hya is 
a particularly fine example, for the higher mass star lies redward of the TAMS in models 
with incomplete mixing physics, a situation which metallicity cannot help. More mixing is 
necessary for good agreement. We do not optimize our models by varying composition. It 
is changed only when an abundance estimate appears in the literature. 

Further insight into potential pitfalls can be found in the binary sample. An increase in 
metallicity moves the tracks in the opposite sense of more complete mixing. It is possible to 
achieve an equally good fit with low metallicity and incomplete mixing or higher metallicity 
and more realistic physics. The practice of making stellar abundance determinations by 
fitting evolutionary tracks is dangerous unless the physics in the code is very well tested 
independently. It is vital that accurate spectroscopic abundance determinations be made for 
stars used as test cases of stellar evolution, particularly eclipsing binaries. 

The pre-MS systems identify an area where our physics is still inadequate. Our predic- 
tive accuracy for these systems is not satisfactory. We must make an ad hoc adjustment to 
the mixing length in order to get large enough model radii. This tells us that our description 
of convection is insufficiently physical. Further multidimensional simulations of envelope 
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convection in low surface gravity stars is necessary to resolve this problem. 

When coupled with the observational tests of light element depletion and turnoff ages in 
young clusters in Young et al. (2003) we explore the performance of TYCHO on stars with 
both convective and radiative cores and convective envelopes of various sizes on the pre-MS 
and main sequence. All of these tests are performed with the same physics. No changes are 
made to the mixing or composition in order to improve our agreement with the observations. 
We find a strong increase in the predictive (as opposed to calibrated) accuracy throughout 
this range of conditions. This is of course a small sub-set of problems in stellar astrophysics, 
but wc have increasing confidence in extending the approach to problems in other areas of 
stellar evolution. In the future we plan to examine the impact of this new approach on 
AGB nucleosynthesis, nucleosynthesis in very low metallicity and evolved massive stars, and 
the evolution of extremely massive stars which become luminous blue variables and SNIb/c 
progenitors. 

This work was supported in part by the DOE (DE-FG03-98DP00214/A001 and the 
ASCII FLASH center at the University of Chicago). We wish to thank J. Liebert, E. Ma- 
majek, C. Meakin, and C. Charbonnel for helpful discussions. We would especially hke to 
thank our scientific editor for his efforts in getting this paper refereed and published in a 
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Table 3. Observed parameters for selected binary systems.'^ 



System 


P(d) 


Star 


Spcct. 


Mass/M0 


Radius/_R0 


logg(cm/s^) 


logre(K) 


log L/Lq 


EM Car 


3.41 


A 


08V 


22.3 ± 0.3*' 


9.34 ± 0.17 


3.864 ± 0.017'^ 


4.531 ± 0.026 


5.02 ± 0.10 


HD97484 




B 


08V 


20.3 ± 0.3*" 


8.33 ± 0.14 


3.905 ± 0.016'" 


4.531 ± 0.026 


4.92 ± 0.10 


V478 Cyg 


2.88 


A 


09.5V 


16.67 ±0.45 


7.423 ± 0.079 


3.919 ± 0.015 


4.484 ± 0.015 


4.63 ± 0.06 


HD193611 




B 


09.5V 


16.31 ± 0.35 


7.423 ± 0.079 


3.909 ± 0.013 


4.485 ± 0.015 


4.63 ± 0.06 


CW Cep 


2.73 


A 


B0.5V 


12.9 ±0.1"^ 


5.685 ±0.130 


4.039 ± 0.024"^ 


4.449 ± 0.011"^ 


4.26 ± 0.06"^ 


HD2 18066 




B 


B0.5V 


11.9 ±0.1"= 


5.177 ±0.129 


4.086 ± 0.024"= 


4.439 ± 0.011"^ 


4.14 ±0.07"" 


QX Car 


4.48 


A 


B2V 


9.267 ±0.122 


4.289 ± 0.091 


4.140 ± 0.020 


4.395 ± 0.009"* 


3.80 ± 0.04"= 


HD86118 




B 


B2V 


8.480 ± 0.122 


4.051 ± 0.091 


4.151 ± 0.021 


4.376 ± 0.010"* 


3.67 ±0.04*= 


CV Vel 


6.89 


A 


B2.5V 


6.100 ± 0.044 


4.087 ± 0.036 


4.000 ± 0.008 


4.254 ± 0.012"* 


3.19 ± 0.05 


HD77464 




B 


B2.5V 


5.996 ± 0.035 


3.948 ± 0.036 


4.023 ± 0.008 


4.251 ± 0.012"^ 


3.15 ± 0.05 


U Oph 


1.68 


A 


B5V 


5.198 ± 0.113 


3.438 ± 0.044 


4.081 ± 0.015 


4.211 ± 0.015"* 


2.87 ±0.08"= 


HD156247 




B 


B6V 


4.683 ± 0.090 


3.005 ± 0.055 


4.153 ±0.018 


4.188 ±0.015"* 


2.66 ± 0.08"= 


C Phe 


1.67 


A 


B6V 


3.930 ± 0.045 


2.851 ± 0.015 


4.122 ± 0.009 


4.149 ± 0.010"* 


2.46 ± 0.04"= 


HD6882 




B 


B8V 


2.551 ± 0.026 


1.853 ± 0.023 


4.309 ± 0.012 


4.072 ± 0.007"* 


1.78 ± 0.04'^ 


IQ Per 


1.74 


A 


B8V 


3.521 ± 0.067 


2.446 ± 0.026 


4.208 ± 0.019 


4.111 ± 0.008"* 


2.17 ±0.03"= 


HD24909 




B 


A6V 


1.737 ±0.031 


1.503 ±0.017 


4.323 ± 0.013 


3.906 ± 0.008"* 


0.93 ± 0.04"= 


PV Cas 


1.75 


A 


B9.5V 


2.815 ± 0.050'* 


2.297 ± O.OSS'* 


4.165 ± o.oie"* 


4.032 ± 0.010"* 


1.80 ± 0.04"= 


HD240208 




B 


B9.5V 


2.756 ± 0.054d 


2.257 ± 0.035'* 


4.171 ±0.016<* 


4.027 ±0.010"* 


1.77 ±0.04"= 


AI Hya 


8.29 


A 


F2m 


2.145 ± 0.038 


3.914 ±0.031 


3.584 ±0.011 


3.851 ± 0.009"* 


1.54 ±0.02"= 


+0°2259 




B 


FOV 


1.978 ± 0.036 


2.766 ± 0.017 


3.850 ±0.010 


3.869 ± 0.009"* 


1.31 ± 0.02"= 


VV Pyx 


4.60 


A 


AlV 


2.101 ± 0.022 


2.167 ±0.020 


4.089 ± 0.009 


3.979 ± 0.009"* 


1.54 ±0.04 


HD71581 




B 


AlV 


2.099 ± 0.019 


2.167 ± 0.020 


4.088 ± 0.009 


3.979 ± 0.009"* 


1.54 ±0.04 


RS Cha 


1.67 


A 


A8V 


1.858 ± 0.016 


2.137 ± 0.055 


4.047 ± 0.023 


3.883 ± 0.010"* 


1.14 ± 0.05"= 


HD75747 




B 


A8V 


1.821 ±0.018 


2.338 ± 0.055 


3.961 ±0.021 


3.859 ±0.010"* 


1.13 ±0.05"= 


EK Cep 


4.43 


A 


A1.5V 


2.029 ± 0.023 


1.579 ±0.007 


4.349 ± 0.010 


3.954 ±0.010 


1.17 ±0.04 


HD206821 




B 


G5Vp 


1.124 ±0.012 


1.320 ±0.015 


4.25 ±0.010 


3.756 ± 0.015 


0.19 ±0.07 


TZ For 


75.67 


A 


F6IV 


1.949 ±0.022 


3.96 ±0.09 


3.532 ±0.02 


3.803 ± 0.007 


1.36 ± 0.06 


HD20301 




B 


G8III 


2.05 ±0.06 


8.32 ±0.12 


2.91 ± 0.017 


3.699 ± 0.009 


1.59 ±0.04 


MY Cyg 


4.01 


A 


FOm 


1.811 ±0.030 


2.193 ± 0.050 


4.007 ±0.021 


3.850 ±0.010"* 


1.03 ±0.04= 


HD193637 




B 


FOm 


1.786 ±0.025 


2.193 ±0.050 


4.014 ±0.021 


3.846 ± 0.010"* 


1.02 ± 0.04"= 


PV Pup 


1.66 


A 


A8V 


1.565 ±0.011 


1.542 ±0.018 


4.257 ±0.010 


3.870 ± O.OIS 


0.81 ± 0.08"= 


HD62863 




B 


A8V 


1.554 ±0.013 


1.499 ±0.018 


4.278 ±0.011 


3.870 ± O.OIS 


0.79 ± 0.08"= 


DM Virf 


4.67 


A 


F7V 


1.454 ±0.008 


1.763 ± 0.017 


4.108 ±0.009 


3.813 ± 0.007 


0.67 ±0.03*= 


HD123423f 




B 


F7V 


1.448 ±0.008 


1.763 ± 0.017 


4.106 ±0.009 


3.813 ±0.020 


0.67 ±0.03*= 


V1143 Cyg 


7.64 


A 


F5V 


1.391 ±0.016 


1.346 ± 0.023 


4.323 ±0.016 


3.820 ± 0.007"* 


0.49 ± 0.03"= 


HD185912 




B 


F5V 


1.347 ±0.013 


1.323 ±0.023 


4.324 ±0.016 


3.816 ± 0.007"* 


0.46 ± 0.03"= 


UX Men 


4.18 


A 


F8V 


1.238 ± 0.006 


1.347 ±0.013 


4.272 ± 0.009 


3.789 ± 0.007S 


0.38 ± 0.03"= 


HD37513 




B 


F8V 


1.198 ± 0.007 


1.274 ±0.013 


4.306 ± 0.009 


3.781 ± 0.007S 


0.32 ± 0.03"= 



''Detailed references and discussion may be found in (Andersen 1991). 
''Stickland, Lloyd, & Corcoran (1994). 



"^Stickland, Koch & Pfeiflfer (1992). 
<^Ribas et al. (2000). 

■^Adjusted here for new T^ff and R determinations. 
'Latham et al. (1996). 
KLastennet & Valls-Gabaud (2002). 
^Hillenbrand & White (2004). 
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Table 4. Results for selected binary systems. 



System 


Star 


Mass log R/Rq 


logTe 


log L log Age (yr) 


z 




hiM Oar 


A 

A 




U. 


.yoy 


4. 


,OzD 


A 

4, 


OOfi 

.yyo 


D.U ( 


0.0189 


0.16 


TTT~\r\T A a A 

HDy7484 


B 


20.30 


0, 


,926 


4, 


,523 


4, 


.898 


6.076 






V478 Cyg 


A 


16.71 


0, 


,881 


4, 


,476 


4, 


,618 


6.310 


0.0189 


0.93 


HD193611 


B 


1 fi Q1 
iD.Ol 


U. 


,ODO 


4, 


A"! A 
,4/4 


A 

4, 


.OoU 


ft Ql O 
D.olz 






CW Cep 


A 


1 o nn 

iz.yu 


U. 




A 

4. 


A'i'7 

,4o i 


A 

4, 


OOQ 

.zzo 


a AC\A 
D.4U4 


0.0189 


0.98 


rlUzloUOD 


B 


1 1 on 
11. yu 


U. 


700 


A 

4. 


AO A 
,4Z4 


A 

4, 


OOQ 


D.4UO 






QX Car 


A 




U. 


,D4U 


4, 


QTO 


Q 

o 


, / Zl 


D.Ool 


0.0189 


0.30 




B 


8.480 


0, 


,602 


4, 


,354 


3, 


,576 


6.563 






CV Vcl 


A 


6.100 


0, 


,609 


4, 


,255 


3, 


,193 


7.295 


0.0189 


0.08 


HD77464 


B 


o.yyo 


U. 


.OUz 


4. 


,z01 


Q 

o 


1 c:o 

.loy 


'7 OOO 

( .zyy 






U Oph 


A 


K 1 OQ 


U. 


.ooo 


4. 


001 

,zzl 


o 

z. 


OOfi 


'7 Q'TO 

i .6 * y 


0.0189 


0.03 


HD15d247 


B 


A fiQQ 


U. 


AQA 


A 

4. 


,iy i 


n 

z, 


700 


'7 QQO 

i .oo\j 






C, Phc 


A 


Q fiQn 

o.you 


U. 


/I c; 
,400 


4. 


,lo i 


z. 


A OO 

.4yu 


1 VOQ 


0.013 


4.03 


rIDoooz 


B 


2.551 


0, 


,261 


4, 


,055 


1, 


,693 


7.728 






IQ Per 


A 


3.521 


0, 


,382 


4, 


,124 


2, 


,211 


7.547 


0.0189 


3.07 


HD24909 


B 


1. / o / 


U. 


. loU 


Q 
O. 


oi Q 

,yio 


n 
U. 


.yoo 


1 ^AT 
( .04/ 






PV Cas 


A 


O QO'7 


U. 


.o4;5 


4. 


OQ? 


1, 


'70'7 

. ( y ( 


ft /loo 


0.0189 


1.00 


rlU24uzuo 


B 




U. 


.OO ( 


A 

4. 


007 

,Uz ( 


1, 


. M Z 


fi Am 

D.4yi 






A I Hya 


A 


O 1 /I 
Z. 140 


U. 


,DUZ 




,Ov54 


1 , 


/I OO 


o.OOO 


0.023 


4.12 


+0 2259 


B 


1.978 


0, 


.434 


3, 


,877 


1, 


.329 


8.558 






V V Pyx 


A 

A 


2.101 


0. 


.339 


3. 


,981 


1, 


.555 


8.610 


0.007 


0.33 


HD71581 


B 




U. 


QQO 

.ooy 


Q 
O. 


,y<5u 


1, 


.OOO 


O.OIZ 






RS Cha 


A 


1.858 


0, 


.317 


3. 


,903 


1, 


.198 


6.866 


0.0189 


2.78 


ilD7o747 


B 


1.821 


0, 


.358 


3. 


,880 


1, 


.189 


6.867 






EK Ccp 


A 


2.029 


0, 


,209 


3, 


,968 


1, 


.242 


7.357 


0.0189 


5.13 


HD206821 


B 


1.124 


0, 


.108 


3, 


,749 


0, 


.165 


7.357 






TZ For 


A 


1.95 


0. 


.552 


3. 


,818 


1, 


.428 


8.713 


0.024 


0.77 


HD20301 


B 


2.05 


0, 


.920 


3, 


,700 


1, 


.331 


8.724 






MY Cyg 


A 


1.811 


0, 


.337 


3. 


,867 


1, 


.095 


8.649 


0.023 


3.79 


HD193637 


B 


1.786 


0, 


.327 


3. 


,865 


1, 


.066 


9.651 






PV Pup 


A 


1.565 


0, 


,184 


3, 


,872 


0, 


,810 


8.520 


0.0189 


0.23 


HD62863 


B 


1.554 


0, 


,180 


3, 


,871 


0, 


.795 


8.527 






DM Vir 


A 


1.460 


0, 


,249 


3. 


,816 





.714 


9.149 


0.023 


0.41 


HD123423 


B 


1.454 


0. 


,243 


3. 


,817 


0, 


.706 


9.149 






V1143 Cyg 


A 


1.391 


0. 


,128 


3. 


,826 


0, 


.515 


8.739 


0.0189 


1.48 


HD185912 


B 


1.347 


0, 


,109 


3. 


,819 





.446 


8.754 






UX Men 


A 


1.238 


0, 


,134 


3, 


,795 


0, 


.400 


9.542 


0.021 


1.56 


HD37513 


B 


1.198 


0, 


,096 


3, 


,795 


0, 


.323 


9.544 
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Table 5. Apsidal comparisons for selected binary systems. 

System Star Mass -logfc, (fcai?-') " P/UcL P/Ugr P/Ucl+gr P/Uqbs 

EM Car A 22.3 1.920 347.0 2.091 0.275 2.37 2.2 ±0.3 

B 20.3 1.935 267.5 

V478 Cyg A 16.67 1.935 160.6 2.97 0.223 3.19 3.0 ± 0.3 

B 16.31 1.917 140.0 

CWCep A 12.90 1.878 49.50 1.46 0.178 1.63 1.640 ±0.014 

B 11.90 1.894 32.11 

QXCar A 9.267 1.898 13.97 0.156 0.170 0.326 0.340 ± 0.006 

B 8.480 1.936 10.58 

U Oph A 5.198 2.124 2.687 1.85 0.0827 1.93 2.2 ± 0.3 

B 4.683 2.111 1.571 

C Phe A 3.930 2.308 1.025 0.894 0.0624 0.956 1.03 ±0.15 

B 2.551 2.333 0.2018 

IQ Per A 3.521 2.278 0.4478 0.335 0.0553 0.410 0.40 ± 0.03 

B 1.737 2.416 0.0401 

PVCas A 2.815 2.149 0.2012 0.414 0.0572 0.499 0.510 ±0.011 

B 2.756 2.221 0.2375 



''Radii R in solaj: units. 
''Multiply tabular value by 10"". 



